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Abstract 

We study the electronic band structure for a model one-dimensional periodic 
potential in the presence of a spacially homogeneous laser field. The statistical 
properties of the energy bands depend on the coupling between the crystal 
and the laser field, going from Poisson to Wigner-Dyson (GOE) and back to 
Poisson as the coupling increases. We argue that the classical dynamics of 
this system resembles that of a periodically driven pendulum. We find that 
the chaotic regime is not restricted solely to high-lying bands and should thus 
be of easier access to optical experiments. 
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In the past decade, many model systems exhibiting, what conventionally became known as 
quantum chaos have been discussed in the literature.chi Among those, the kicked rotatorS^B 
and the Hydrogen atom in the presence of a uniform magnetic fielda have been studied in 
great detail, both theoretically and experimentally. Formally, one speaks of quantum chaotic 
behavior when the system dynamics in the classical limit is very irregular or completely 
chaotic. For a quantum system whose classical dynamics is regular one usually finds that 
energy levels are uncorrelated and therefore obey a Poissonian statistics. On the other hand, 
if the system is fully chaotic in the classical limit, its energy levels tend to repel each other 
strongly. In this case the statistical properties of the spectrum coincide with those drawn 
from the Gaussian ensembles of random matrix theory (RMT).I 

Recently, Mucciolo and coworkers showedS that the electronic spectra of crystalline 
solids (e.g.. Si) also show the universal signatures of quantum chaotic behavior. Their 
analysis was restricted to high-lying bands, well above the Fermi energy. Since crystals 
are fixed by nature, no external parameter exists to modify the system and therefore try to 
enhance the irregularity in its classical dynamics. This is in contrast with the two previously 
mentioned systems, where chaos is controlled by an external parameter, as for example 
the strength of a time-periodic impulse or the intensity of a magnetic field. 

In this work we present a new model exhibiting spectral signatures of chaos based on 
electrons in a periodic potential. In our case a monochromatic laser field is the main respon- 
sible for the irregular electron dynamics. The intensity of the laser field can be regarded 
as an external parameter changing the behavior of the system. We believe that this mech- 
anism could increase the chances of experimentally observing chaos in crystals by bringing 
universal statistics to low-lying bands as well. By looking at the energy level statistics of 
a very simple 1-D periodic potential we were able to identify three regimes. These suggest 
that the system dynamics moves from regular to chaotic and then back to regular as the 
field intensity is increased. Our formalism is based on the dressed band approachjll which 
we briefly review below. 

According to Bloch's theorem, the wave function for an electron in a 1-D periodic poten- 
tial can be labelled by a momentum k, —n/d < k < ir/d, where d is the lattice parameter. 
This allows the Hamiltonian to be decoupled into well-defined k (reduced) components. In- 
clusion of the interaction with a monochromatic laser field leads to the following reduced 
Hamiltonian: 

H = Hk + + Hint , (1) 
where, in atomic units {h = e = m = 1), 
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H,r^t = aAp+ —A^ . (4) 

In the above equations, V{z) denotes the lattice potential, (a) is the creation (annihila- 
tion) field operator, a = 1/137 is the fine structure constant, and A = ^^J2tt/ujQ (a + a"*"), 
where Q is the quantization volume. The periodicity of the lattice potential allows it to be 
written in terms of a reciprocal lattice sum over = 2'n'i/d, as indicated in Eq. (^. Here 
the field polarization is taken along the z-direction. 

By performing adequate unitary transformations in terms of a basis of eigenstates \G, n) 
of the momentum p and the field Hamiltonian if-y, H can be transformed into a Floquet 
Hamiltonian Hp- One can then show0 that the spectrum of H follows from the diagonal- 
ization of the Bloch-Floquet matrix 



{G\n'\HF- NQUj\G,n) 



{G + kf 
nuj + -^^ — ^ 



^G,G'^n,n' + Jn'-n(y—^{G - G')^ Vg-G' , (5) 

where Aq = ujVL is the amplitude of the classical vector potential and J„ is the 

Bessel function of order n. It should be noted that n in the state |G,n) denotes an integer 
(positive or negative) that describes the deviation from the average number of photons Nq, 
which is supposed to be large. Note that the matrix defined by Eq. (|^) is time-reversal 
symmetric. 

We present results for a model potential where V{Gii) = —a for —3 < ^ < 3 and zero for 
other i values. This is a truncated Kronig- Penney attractive model. We take the energy unit 
to be e = Gl/2 and the potential strength a = 0.25e. The laser intensity is characterized by 
the dimensionless parameter x = aAoGi/u. Equation (H) allows for a simple interpretation 
of the spectrum in this problem: In the absence of the laser field, the electronic energy levels 
consist essentially of the free-electron parabola folded into the first Brillouin zone (BZ), plus 
gaps opening around the crossing points {k = and ±Gi/2). The gap amplitude and the 
distortion of the bands with respect to the free-electron parabola depend on the potential 
strength a. Inclusion of the field produces a dressing effect in the bands, which may be 
described by replica-bands translated by an amount nu. The interaction Hint causes an 
anti-crossing whenever two noninteracting replica-bands cross. Clearly, this occurs for all k 
in the BZ; for a fixed value of u, the gaps widen as the intensity x increases. 

In Fig. we show the spectrum obtained for —10 < n < 10, u = 0.3e, and field intensity 
X = 10. We avoid regions at the extremes of the energy spectrum, where only replicas of 
the lowest or higher bands are present (due to the truncation in n adopted in our numerical 
treatment) and limit our analysis to the energy region 1 < E/e < 5, where a large number 
of anti-crossings occur. For the parameters chosen here, this energy region is representative 
of the complete basis, — oo < n < oo. The laser frequency u is chosen to maximize the 
mixing between the dressed bands. Note that the spectrum in Fig. ^ reveals an intricate 
level repulsion structure which resembles those characteristic of quantum chaotic regimes in 
other systems.ii'ill2l 

A quantitative identification of quantum chaos consists primarily in unfolding the energy 
spectrum (to attain a constant density of states) and then analyzing the short and long-range 
statistical correlations present. For spinless, time-reversal symmetric systems with a complex 
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dynamics, predicts that the distribution of nearest-neighbor level spacings (NNS) 

should follow closely the Wigner surmise for the Gaussian orthogonal ensemble (GOE), 

pGOE^^^ ^ g-..V4 ^ (6) 

where the level spacing s is given in units of the mean level spacing. On the other hand, the 
energy levels of most integrable systems are essentially uncorrelated (one important excep- 
tion is the harmonic oscillator), in which case the NNS distribution is a simple exponential, 
namely, 

pPoisson^^^ = e-' . (7) 

Long-range correlations can be diagnosed through the least square deviation A^fL) of the 
number of levels in a given energy interval of length L from a constant slope." When L 
is given in units of the mean level spacing, A3(L) is expected to fall between two extreme 
limits: linear (L/15) for uncorrelated levels and constant (1/12) for fully rigid, equally spaced 
spectra. The GOE prediction! is an intermediate function which, for large L, behaves as 
AGOE(^) ~ (1/7^2) ln(27rL) - 0.007. 

In Fig. 1^ and || we present our results for the NNS distributions and A3, respectively, 
obtained for increasing values of the field intensity parameter x. The data correspond to a 
sampling over 100 equidistant k values along < /c/Gi < 0.5. The level density remained 
essentially constant on the energy interval considered and no unfolding was carried out. 
The resulting statistics involved over 5,200 calculated levels for each x. Let us concentrate 
initially on Fig. ^. For a; = the uncorrelated character of the energy levels and the 
integrability of the Hamiltonian in Eq. (|1]) is reflected in the agreement between the data 
and the Poisson law. As x increases, a weak level repulsion appears. Note that x = \ 
corresponds to a crossover regime in which neither of the distributions in Eqs. (0) and (^ 
agrees with the data. However, for x = 10 the Wigner surmise provides a very accurate 
fit, suggesting a chaotic regime at this intensity. Further increase in x causes the data to 
change, deviating from the Wigner surmise. For x = 1000, the spectrum tends again to a 
Poisson statistics, which reflects the exact integrability in the limit of infinite field intensity 
as well. It is important to note that the NNS is not a delta-function centered at uj when 
X ^ 00: The periodic potential, although relatively weak, folds the sequence of equidistant 
parabolas back onto the first BZ, causing the NNS to be Poissonian. 

A similar interpretation goes for Fig. ^. The maximal spectral rigidity occurs around 
X = 10. In this figure, however, we see a marked deviation from either the Poisson (x = 
and 1000) or the GOE [x = 10) predictions beyond intervals of length Lmax ~ 4. From 
Fig we may notice a tendency of levels to form clusters of four. This effect is well known 
in the literature:t3'ta In a semiclassical language, it indicates the breakdown of the universal 
regime and the appearance of system-dependent contributions, which are related to very 
short periodic orbits. We have checked that //max is always approximately equal to the 
laser frequency u (in units of the mean level spacing). The periodic orbit associated to this 
frequency is the oscillation of the laser field in its phase space.0 The reason why this orbit is 
always the shortest possible one and its signature survives different values of x is related to 
the nature of our model. For the Kronig- Penney model there is only one bound state {E < 0) 
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in the spectrum. As a result, in the energy range where the A3 was obtained, the electron 
classical motion can only be bounded (and therefore periodic) through the interaction with 
the laser field. This always leads to a periodicity T > li^juj and the intrinsic period of 
oscillation of the laser field remains the shortest one. 

As noted before, for a 3-D periodic potential the electronic spectrum at low-symmetry 
regions of the BZ manifests quantum chaotic behavior. In the equivalent 1-D problem, no 
chaotic behavior is expected, since the electrons have only one degree of freedom iz) and 
the problem is classically integrable. The introduction of a laser field adds another degree of 
freedom to the system and makes the problem nonintegrable. As indicated by our numerical 
results, for high enough intensities of the laser field, it is plausible that this perturbation 
causes the electron motion to be strongly irregular, leading to chaos. 

Although we treat the laser field in a second quantized form, our results are equivalent to 
those obtained from a semiclassical, time- dependent description where the field is fixed and 
can be written in the form cos(ti;t). Both descriptions yield the same spectrum (band 
structure) of energies and quasi-energies in the quantized and semi-classical descriptions, 
respectively. The only difference is that, in the semiclassical case, one formally incorporates 
the time-periodic vector potential to the crystal potential, resulting in a space-time (two- 
dimensional) BZ, but with matrix elements for the effective Hamiltonian similar to Eq. (^. 

It is interesting to compare the problem of a laser- illuminated 1-D crystal to the kicked 
rotator .i'S The dynamics in the latter is generated by a Hamiltonian in which the free rotator 
is subject to a time-periodic force whose amplitude is a function of the rotation angle (27r- 
periodic in ff). In the classical counterpart of our model (and in the kicked rotator), an 
external field, space and time periodic, acts on a particle moving in one dimension. The 
difference between these two systems is that, in our case, space and time components are 
decoupled, allowing one of the potentials to be made negligible with respect to the other by 
varying the parameter x. Thus, the x = limit corresponds to an electron in a 1-D crystal 
problem, and a; — > 00 corresponds to a free electron in the presence of the laser field. In 
both limits the problem is exactly integrable. 

No such decoupling is possible in the kicked rotator. Our problem would be more closely 
related to a periodically kicked "pendulum^ in which case the kick amplitude A would be a 
control parameter similar to the laser intensity, and this system should present quantum 
chaos for intermediate values of A. For A = and for A 00 the problem would be 
dominated by the gravitational field or by the kick force respectively, becoming integrable 
and therefore non-chaotic. 

To our knowledge, no other system discussed in the literature has been led from regularity 
into chaos and than back to regularity by increasing a single control parameter. One might 
think that a similar situation would occur for the H atom in a magnetic field since the 
problem is exactly integrable in the limits of zero magnetic field and zero Coulomb potential. 
However, even in the limit of extremely strong magnetic field, the Coulomb interaction can 
not be neglected. In this limit the classical motion of the electron is strongly constrained 
in the field direction and the electronic distribution is no longer spherically symmetric. The 
electronic motion may be approximated by a 1-D truncated Coulomb interaction with the 
nucleus, which has a completely different spectrum from the Landau problem.B Other simple 
model-Hamiltonians, such as the the spin-boson model,lli may display a variety of behaviors, 
but they require variations of several control parameters. 
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Finally, we discuss the experimental implications of our model. Contrary to the 3-D band 
structure problemji the chaotic behavior presented here is not limited to highly excited 
bands. Any region of the spectrum of the 1-D system displays the universal features under 
a laser field of adequate intensity. The same effect should be present in 3-D systems, even 
for bands close to the Fermi energy. The limitation on the long-range universal correlations 
imposed by short periodic orbits should become less stringent in higher dimensions due to 
the increase in the density of states. We believe that optical pump-probe techniques should 
be adequate to investigate the dressing effects presented here. They should reveal universal 
signatures in the frequency and momentum dependence of the dielectric function and the 
optical conductivity, as predicted by Taniguchi and Altshuler.0 
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FIGURES 




FIG. 1. The electronic spectrum in a 1-D periodic potential under a laser field of intensity 
X = 10. Energy bands are shown as a function of the Bloch momentum k. 
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FIG. 2. Histograms giving the nearest-neighbor spacing distributions for the energy spectra of 
Eq. (H) under increasing field intensity (from top to bottom). Sohd and dashed hues correspond to 
GOE and Poisson statistics respectively. Note that, as the field intensity increases, the distribution 
evolves from Poisson-like into GOE-like, and then back to Poisson-like. 
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FIG. 3. Least square deviation statistics A3(L) for the indicated values of the field intensity 
parameter x. The sohd and dashed lines are the GOE and Poisson regimes, respectively. 
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